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I  INTRODUCTION 


The  propagation  of  a  high  current  relativistic  electron  beam  through  a 
gaseous  medium  is  strongly  influenced  by  the  conductivity  of  the  medium.  For 
an  initial  beam  pulse  into  air,  the  conductivity  is  generated  by  collisional 
ionization  of  the  air  by  primary  beam  electrons  and  by  the  numerically  domi¬ 
nant  cascade  or  secondary  electrons;  if  the  beam-generated  inductive  electric 
fields  are  strong  enough,  additional  ionization  is  created  by  the  field-driven 
secondary  electrons.  For  typical  beam  parameters  into  atmospheric  density 
(1  amagat)  air,  electron  densities  in  the  beam  channel  are  on  the  order  of 
10^/cm^,  which  corresponds  to  a  conductivity  (assuming  T„  >  300  K)  of  order 
a  *  10^/s  (units  of  equal  conductivity  are  9  x  10^  s~^  *  1  mho-m-*  = 

10'2  mho-cm-'- ) . 

When  multiple  beam  pulses  are  fired  in  a  sequence,  the  residual  conducti¬ 
vity  in  the  afterglow  air  channel  from  previous  pulses  can  become  large  enough 

to  influence  the  propagation  of  a  subsequent  pulse.  The  critical  value  is  of 
9  —1 

order  o  *  3  x  10  s  ,  which  corresponds  to  room  temperature  electron 
densities  nfi  =  10  /cm  .  Beams  of  interest  will  be  approximately  1  cm  in 
diameter. 

The  objective  of  the  work  described  here  is  to  develop  and  demonstrate  a 

diagnostic  technique  capable  of  measuring  conductivity  levels  of 
-9  -1 

a  <  3  x  10  s  in  a  1-cm-diameter  plasma.  Since  physical  probes  would  most 
likely  perturb  the  beam,  and  would  be  subject  to  damage  as  well,  the  techni¬ 
ques  considered  must  be  nonintrusive ,  which  means  that  they  will  involve  the 
interaction  of  the  plasma  with  a  diagnostic  eletromagnetic  wave  or  field. 

We  first  describe  that  interaction  for  EM  waves  propagated  through  the 
plasma,  and  we  show  sensitivity  limits  for  this  approach  (i.e.,  classical 
interferometry  and  absorptivity  in  the  infrared  through  microwave  frequency 
ranges).  After  showing  that  such  single-pass  measurement  techniques  are 
Inadequate,  we  then  consider  a  multiple-pass  or  standing  wave  technique, 
namely,  a  microwave  cavity  resonance  perturbation  technique,  which  is  well 
suited  for  these  measurements.  Optimization  of  the  cavity  design  and  antici¬ 
pated  performance  characteristics  of  the  system  are  then  described. 
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The  resulting  measurement  system  has  been  used  to  perform  several  pre¬ 
liminary  checkout  measurements  of  the  conductivity  in  the  afterglow  of  air 
excited  by  a  Febetron  706  electron  beam.  These  measurements  give  first  re¬ 
sults  for  room  temperature  laboratory  and  synthetic  air  over  the  pressure 
range  from  10  to  760  torr.  As  these  measurements  are  refined,  they  will  pro¬ 
vide  benchmarks  for  calculating  the  afterglow  plasma  properties  using  air 
chemistry  codes.  Air  chemistry  is  sufficiently  complex  to  necessitate  such 
comparisons  at  room  temperature,  and  even  more  difficult  at  the  elevated 
temperatures  characteristic  of  heated  channels.  We  plan  to  address  this  pro¬ 
blem  in  the  future  by  making  measurements  in  shock-heated  air. 

When  electron  and  ion  densities  are  anywhere  near  comparable  magnitudes, 
the  plasma  properties  are  dominated  by  the  electrons.  However,  in  a  medium 
where  electron  attachment  rates  are  large  compared  with  electron-positive  ion 
recombination  rates,  negative  ion  densities  can  be  orders  of  magnitude  higher 
than  electron  densities  in  the  late  afterglow.  This  might  be  the  case  for 
room  temperature  air  at  atmospheric  pressure,  and  we  find  anomalies  in  our 
1  amagat  data  that  suggest  a  role  of  negative  ions.  While  this  is  of  interest 
of  itself,  it  also  illustrates  the  possibility  that  the  conductivity  could 
rapidly  change  when  fields  from  a  propagating  pulse  interact  with  the  after¬ 
glow  plasma.  It  may  be  possible  to  simulate  this  effect  in  our  experiments  by 
applying  a  short  pulsed  field,  either  DC  or  microwave,  to  the  afterglow  plasma 
in  the  cavity  and  measuring  the  resulting  charge  in  conductivity. 


II  ELECTROMAGNETIC  PROPERTIES  OF  PLASMAS 


This  subject  is  a  large  and  complicated  one,  which  easily  occupies  entire 
books.  A  particularly  useful  source  relevant  to  diagnostics  of  plasmas  is  the 
book  by  Heald  and  Wharton,  Plasma  Diagnostics  with  Microwaves  [HW65].  Here  we 
will  summarize  the  relevant  equations  required  in  our  analyses. 

Complex  Propagation  Constants 

The  scalar  dielectric  constant  of  a  linear  isotropic  medium,  x,  is  the 
ratio  of  permittivities  of  a  medium  to  vacuum: 

<  =  e/e0  ’  (2.1) 

while  the  scalar  conductivity  of  such  a  medium,  a,  is 

J  =  oE  .  (2.2) 

These  two  parameters  can  be  combined  to  form  a  complex  dielectric  constant,  ic, 
as 


<  =  Kr  ”  iKi  “  K  ~  jo/a,e0  *  (2.3) 

where  j  *  /-T  .  Likewise,  one  can  define  a  complex  conductivity,  a: 

a  h  or  +  ja1  -  a  +  ju>(K  -  l)e0  .  (2.4) 

Comparison  of  the  two  expressions  shows  that  they  are  related  by 

"a  a.  o_ 

*  “  Kr  “  3Ki  *  1  ~  J  - -  (1  +  — L  )  -  j  — L  (2.5) 

eow  eo“  Eo“ 


and 


o*o  +  jo.  ■  je  ui(k  -  1)  ■  x.e  u  +  j(<  -  1  )e  u> 

r  i  o  io  J  r  o 


(2.6) 
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Since  is  related  to  dissipation,  it  is  common  in  microwave  measurements  to 
refer  to  the  "loss  factor"  =  *  or/e0ui.  Further,  the  "loss  tangent"  is 

defined  as 


tan  6 


<  e  u)K 
r  o 


(2.7) 


Propagation  of  an  EM  wave  in  a  plasma  can  conveniently  be  described  in 
terms  of  a  complex  refractive  index,  p,  which  is  related  to  the  complex 
dielectric  constant  by 


u  **  u  -  ix  =  •  (2.8) 

The  fringe  shift,  AS,  is  related  to  the  change  in  the  real  part  of  the  index 
as 


AS  =  (pQ  -  p)/X  cm" 


(2.9) 


while  che  attenuation  of  a  probe  beam  is 


a  *  xu/c  cm 


-1 


(2.10) 


In  terms  of  the  complex  conductivity,  the  refractive  indices  can  be  written  as 


(2.11) 


Although  the  parameters  described  above  are  usually  presented  as  approxi¬ 
mations  in  one  limit  or  another  (D.C.,  u>  «  v;  A.C.,  u  »  v;  v  ■  constant), 
all  of  the  parameters  can  be  given  in  complete  generality  through  the  expres¬ 
sion  for  complex  conductivity  in  terms  of  the  momentum  transfer  collision  fre¬ 
quency  v(v)  and  the  electron  velocity  distribution  function  fQ(v)  as  [HW65] 


a 


2  m 

4ti  ne  r  1 

3  m  £  v  +  ju> 


df  (v) 
o 

dv 


3  . 

v  dv 


(2.12) 
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This  expression  can  be  put  in  more  familiar  form  by  defining  a  new  com¬ 
plex  function,  namely,  a  complex  collision  frequency  v  h  vr  +  jv^.  This 
function  can  be  defined  in  terms  of  the  integral  in  equation  (2.12)  for  arbi¬ 
trary  electron  distribution  functions,  but  when  the  distribution  is 
Maxwellian,  it  becomes  explicitly 

c  j9  2 

_  8  m  «  v  +  jto  4  -mv 

V  1  VT  W  L  - 2~T  v  exp  t-lRTi  dv  ‘  (2.13) 

ZKl  0  1  +  V  /u) 


This  parameter  can  be  calculated  for  any  gas  (or  gas  mixture)  for  which  v(v) 
is  known.  Then,  by  introducing  in  addition  the  plasma  frequency,  uip,  where 


2 


to 


P 


nee 

meEo 


(2.14) 


the  real  and  imaginary  conductivities  can  be  written  simply  as 

2 


°r  vr  w 


CQU  <0  (0 


U) 


(2.15) 


(2.16) 


Mote  again  that  this  description  of  the  conductivity  is  completely 
general,  applying  to  arbitrary  values  of  to  ,  to,  and  v.  It  is  not  even 
necessary  that  the  electron  distribution  function  be  Maxwellian;  but  it  must 
be  known  if  equation  (2.12)  is  to  be  evaluated.  However,  it  is  cumbersome  to 
re-evaluate  the  integral  in  equation  (2.13)  for  each  combination  of  u  and  v, 
so  it  is  worth  examining  some  special  cases.  First,  when  to  «  v  (the  D.C. 
limit). 


2  3 

to  to 

+  j  ~ 

V  V 


When  to  »  v  (the  A.C.  limit), 


v  +  jto 

,  2 

1  +  v  /to 


(2.17) 
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v  +  ju 


(2.18) 


^  rj 

Vie  will  therefore  define  (veff)jx;»  ^weff^AC*  anc*  (veff  ^DC’  as  ^°^ows; 


(2.19) 


(2.20) 


(2.21) 


The  effective  collision  frequencies  in  equations  (2.19)  and  (2.20)  were  dis¬ 
cussed  by  Heald  and  Wharton  [HW65],  and  have  been  evaluated  for  several  common 
molecular  gases  by  Itikawa  [ It 7 1 , Xt7 3 ) .  We  introduce  the  third  definition  to 
write  the  simplified  forms  of  the  real  and  imaginary  conductivities  in 
Table  1.  These  effective  collision  frequencies  are  only  a  function  of 
electron  temperature.  Itikawa's  results  can  be  used  to  calculate  values  for 
dry  air  (80%  N2,  20%  Oj)  or  50%  humid  air  (add  8.5  torr  1^0),  with  the  results 
shown  in  Figure  1.  The  values  at  the  A.C.  and  D.C.  limits  are  generally  with¬ 
in  50%  of  each  other,  suggesting  (but  not  proving)  that  use  of  an  average 
value  over  the  entire  frequency  range  is  not  a  bad  approximation. 

For  completeness,  we  can  substitute  the  expressions  from  Table  1  into 
equation  (2.11)  to  arrive  at  limits  for  the  real  and  imaginary  refractive 
indices.  These  limits  are  given  in  Table  2. 


Sensitivity  of  Simple  Pass  Measurements 

To  evaluate  the  sensitivity  of  single-pass  probe  experiments,  parti¬ 
cularly  in  the  far-infrared  and  millimeter  wavelength  ranges,  we  have  calcu¬ 
lated  the  indices  of  refraction  using  the  v  =  constant  expressions  of 

Table  2  and  the  (veff)  values  for  Tg  *»  300  K  (0.025  eV)  from  Figure  1.  The 

DC 

resulting  AS  and  a  values  are  plotted  versus  probe  wavelength  in  Figures  2 
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Table  2 


COMPLEX  REFRACTIVE  INDICES  IN  SPECIAL  CASES 


u)  ^  v  (D»C«) 
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FIGURE  1  EFFECTIVE  COLLISION  FREQUENCIES  IN  AIR  AT  1  amagat 
[From  V.  Itikawa,  Pbvs.  Fluids  16,  831  (1973)1 


Fringes/cm 
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FIGURE  2  FRINGE  SHIFT  DUE  TO  ELECTRONS  IN  A  LOSSY  PLASMA 
(1  amagat,  Te  =  300°  K) 


L0 


and  3  for  dry  and  50%  humid  air  at  1  amagat.  In  the  A.C.  limit  (short  wave- 

o 

length),  AS  and  a  increase  with  X  and  X  ,  respectively,  and  one  tends  to 
assume  that  sensitivity  can  be  increased  by  increasing  the  wavelength  indefi¬ 
nitely;  however,  this  is  clearly  not  the  case:  both  AS  and  o  reach  a  maximum 
when  oj  *  v,  and  subsequent  increases  in  X  decrease  the  sensitivity. 

If  one  assumes  sensitivity  limits  of,  say,  a  1%  change  in  measured  inten- 

_2 

sity  for  a  1-cm-diameter  plasma  at  1  amagat,  then  =  10  /cm  and 

AS^n  a  1*6  x  lO'^/cm.  These  limits  correspond  to  ng  =  lO^/cnr*  in  both 
cases.  In  fact,  the  sensitivity  is  not  likely  to  be  this  good,  since  probe 
beam  diameters  will  be  about  the  same  size  as  the  plasma  diamter,  leading  to 
strong  deleterious  refraction  effects  [KW64,  Ves75,  Ver79,  DK75,  GF71, 

Sh61].  Thus,  a  practical  sensitivity  limit  for  a  single  pass  measurement  is 
about  ng  =  10  /cm  at  1  amagat.  Fringe  shift  measurements  can  be  more 
sensitive  at  lower  densities  since  the  turnover  due  to  terms  in  v  occurs  at 
longer  wavelengths,  but  the  attenuation  is  proportional  to  v,  and  so  it 
becomes  less  sensitive.  Note  that  in  the  range  of  maximum  sensitivity,  both  a 
and  AS  are  functions  of  both  a ^  and  or,  equation  (2.11),  so  it  is  necessary  to 
measure  both  a  and  AS  to  determine  either  or  or. 

Multiple  Pass  Measurements 

Greater  sensitivity  can  be  achieved  by  increasing  the  interaction  path- 
length  between  the  probe  beam  and  the  plasma,  in  this  case  by  probing  with 
multiple  passes  through  the  plasma.  At  optical  through  millimeter  wave¬ 
lengths,  this  would  be  accomplished,  for  example,  by  switching  from  a  single 
pass  Mach-Zehnder  interferometer  to  a  multiple  pass  Fabry-Perot  interfero¬ 
meter.  However,  such  an  approach  is  not  entirely  satisfactory,  since  the 
response  of  the  transmission  of  a  Fabry-Perot  interferometer  to  an  absorption 
is  highly  nonlinear,  and  since  refraction  effects  still  outweigh  interference 
or  absorption  effects.  At  microwave  frequencies,  multiple  pass  sensitivity  is 
achieved  by  using  a  cavity  that  completely  contains  the  fields,  thus  elimi¬ 
nating  refraction  effects.  Furthermore,  the  response  of  the  cavity  to  pertur¬ 
bations  is  well  defined.  Me  devote  the  remainder  of  this  report  to  the  micro¬ 
wave  cavity  approach. 


1  1 


Ill  THE  MICROWAVE  RESONANT  CAVITY  TECHNIQUE 


Description 

The  properties  of  microwave  cavities  are  described  in  several  books, 
including  Montomery's  Technique  of  Microwave  Measurements  (Mo47).  Basically, 
a  cavity  supports  transverse  electric  (TE)  and  transverse  magnetic  (TM)  modes, 
each  with  a  large  number  of  resonant  frequencies.  For  a  right  cylindrical 
cavity,  the  modes  are  further  characterized  by  the  subscripts  "£mn,”  where  l 
is  the  number  of  standing  waves  in  the  azimuthal  direction,  ra  is  the  number  of 
half  waves  in  the  radial  direction,  and  n  is  the  number  of  half  waves  in  the 
axial  direction.  The  resonant  wavelengths  are 
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(3.1) 


where  the  are  the  m'th  roots  of  the  Bessel  functions  J^(x).  A  few  of  the 
lowest  mode  frequencies  are  shown  in  Figure  4.  Each  of  these  resonances  has 
an  intrinsic  unloaded  cavity  Q  ■  Qu,  where  Q  is  defined  as  2tt  x  the  ratio  of 
energy  stored  to  the  energy  loss  per  cycle  of  oscillation  of  the  resonant 
standing  wave  (due  to  resistive  losses  of  the  surface  currents  in  the 
cavity).  Again,  for  right  cylindrical  cavities,  Qu  is  related  to  the  skin 
depth  6  by 
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where  p  ■  nx/2.  The  skin  depth  is  further  related  to  the  resistivity  p  as 

xp _ ,  V2. 
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(tO)2  X  10'20  -  (Hz  -  cm) 


A  high  Qu  is  achieved  by  choosing  a  cavity  with  high  volume-to-surface  ratio 
and  by  plating  the  interior  surface  with  a  low  resistivity  material  like 
silver. 

Both  the  resonant  frequencies  and  the  Qs  of  a  cavity  are  perturbed  by 

introducing  a  material  with  a  *  <  ,  where  <_  is  the  dielectric  constant  of 

o  ° 

air.  Slater  [S146]  derived  a  perturbation  formula  as 
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When  the  perturbation  is  small  enough  that  the  electric  fields  retain  their 
initial  values,  E  »  E  ,  the  perturbation  formula  can  be  written 
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(3.7) 


where  we  have  now  identified  the  perturbation  as  a  plasma.  The  last  expres¬ 
sion  shows  clearly  that  if  the  area  of  the  plasma  is  known  so  that  the  ratio 
of  integrals  can  be  evaluated,  then  the  real  and  imaginary  conductivities,  o_ 
and  Oj,  can  be  determined  by  measuring  a(— )  and  Aw/w,  respectively.  This  is 
in  contrast  to  the  single  pass  probe  measurement,  described  in  the  previous 
section,  where  both  AS  and  a  had  to  be  measured  tc-  determine  crr  or  <jj. 


Me  emphasize  again  that  in  the  resonant  cavity  perturbation  technique, 
the  change  in  cavity  Q  is  directly  proportional  to  or,  and,  furthermore,  if 
u>  «  v,  this  gives  the  D.C.  conductivity,  or  *  0d.c.#  71118  18  the  parameter 

of  primary  Interest  in  these  measurements. 


The  Slater  formula  requires  corrections  when  the  perturbations  become 
large.  In  our  case,  our  right  cylindrical  cavity  must  have  holes  in  the  ends 
to  pass  the  electron  beam,  and,  in  the  Febetron  experiments  reported  here,  a 
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quartz  tube  was  inserted  the  length  of  the  cavity  to  confine  the  plasma  to  a 
known  area.  Corrections  for  such  effects  have  been  derived  for  Aw/w  versus 
Oj.  [Th63,  Th65,  Lu68],  but  apparently  not  for  a(qJ  versus  or.  Therefore,  we 
base  our  results  to  date  on  the  Slater  formula,  with  the  plan  to  calibrate  the 
correction  using  liquids  of  known  conductivity  [BF49,  SL461. 

There  are  two  ways  to  determine  the  A(l/Q)  change  due  to  a  plasma  pertur¬ 
bation.  The  first  is  to  measure  the  cavity  transmission  at  several  frequen¬ 
cies  in  the  vicinity  of  the  resonance  and  to  cross  plot  the  power  transmitted 
at  finite  time  intervals  following  the  e-beam  pulse.  Then,  by  definition, 

Q(t)  *  fQ ( t ) / Af ( t ) •  The  second  is  based  on  the  relative  value  of  the 
transmission  itself.  For  the  empty  cavity,  the  losses  giving  rise  to  the 
loaded  are 
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where  0j  and  0^  are  defined  in  terms  of  the  ratios  of  the  coupling  Qs  to  Qu. 
The  power  transmission  through  the  cavity  is  then  [Ha63] 

4  0i  02 

TC  “  2  2  w  -  u)  T  (3.9) 

(1  +  01  +  02 )  +  *QU 

r 

At  the  resonant  frequency  w  *  wr,  this  simplifies  to 


t  0  ■ 
c 


<V°ext> 


(3.10) 


When  a  plasma  is  introduced,  the  unloaded  Q  is  perturbed  by 


(3.11) 
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while  the  external  0  stays  constant 
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xt 
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(3.12) 


It  can  be  shown  that  the  ratio  of  cavity  transmissions  with  and  without  a 
plasma  is  then 


(3.13) 


If  both  measurements  are  done  at  the  resonance  peak,  this  becomes  simply 


(3.14) 


Thus,  a  measurement  of  the  cavity  transmission  at  the  resonance  peak  suffices 
to  give  &(1/Q)(t)  and  therefore  ar(t). 


Cavity  Optimization 

Maximum  sensitivity  to  a  centerline  plasma  is  achieved  by  choosing  a 
™CSMO  m0^e»  f°r  which  the  fields  are 


Er  *  Eg  *  0 


(3.15) 


and 


EZ  *  Jo  <X*mr/Rc>  •  (3.16) 

The  overlap  function,  which  is  the  ratio  of  integrals  in  equation  (3.7),  then 
has  the  values  shown  in  Figure  5.  While  this  figure  suggests  that  maximum 
sensitivity  is  achieved  by  using  a  small  radius  cavity  and  a  high  order  mode 
(m  >  3),  a  more  detailed  look  at  cavity  optimization  is  appropriate. 
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We  begin  the  analysis  with  the  equation  for  cavity  transmission,  equation 
(3.10).  Taking  the  partial  derivative,  and  using  equation  (3.11),  yields 

3T  . 

-f - 2V(q)  ■  -2V(o)  <3-17> 

C 

Optimization  of  the  cavity  now  consists  of  maximizing  the  product  of  Qu  and 
A(l/Q).  For  a  right  cylindrical  cavity  with  TM^mg  modes,  Qu  was  given  in 
equation  (3.3).  Since  the  skin  depth,  5,  is  proportional  to  A1'  ,  6  =  6qa  , 

the  complete  expression  can  be  written 


ou  = _ _ 

U  5  /X  it(2  +  D/L) 
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(3.18) 


To  evaluate  &(1/Q),  we  find  that  the  overlap  function  can  be  approximated  by 
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Combining  equations  (3.18)  and  (3.20)  into  (3.17)  with  the  additional  note 
that  xjj,m  ”  2ttRc/A  yields  the  desired  result 
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We  want  a  cavity  that  has  the  maximum  3TC/TC  for  a  given  or  and  Rp.  We 
see  that  the  sensitivity  is  independent  of  the  mode  (within  the  TM^q  manifold) 
and  the  cavity  diameter,  and  depends  only  on  oQ,  D/L,  and  A^^  (obviously,  a 
and  mode  determine  the  diameter).  5Q  and  D/L  should  be  minimized  to  achieve 
maximum  sensitivity;  this  maximizes  Qu,  which  is  desirable  as  can  readily  be 
seen  from  equation  (3.17). 

We  note  that  the  wavelength  dependence  suggests  as  short  a  wavelength  as 
possible.  However,  there  are  some  limitations  on  reducing  that  parameter. 
Guidelines  in  the  literature  suggest  that  the  length  of  the  cavity  be  at  least 
4  times  the  hole  diameter  [SL46,  Th63,  Th65]  and  that  the  cavity  radius  be 
about  10  times  the  plasma  diameter  [Lu68],  Furthermore,  A  should  be  greater 
than  ,IDaperCure  to  avoid  radiation  losses  out  the  holes.  Finally,  D/L  >  1 
facilitates  operation  on  the  TMQfflQ  modes  without  interference  from  higher 
order  modes. 

We  have  chosen  t>aperture  =  3  cm  to  accomodate  electron  beam  passage 
without  interference  on  ETA  and  ATA  experiments.  We  have  compromised  on  the 
above  conflicting  requirements  by  choosing  L  =  11.94  cm  and  D  =  14.73  cm 
operating  in  the  TMq^q  mode.  The  calculated  resonance  frequency  in  this  case 
is  1.558  GHz,  while  the  measured  value  is  1.5624  GHz  (A  ■  19.2  cm).  The 
cavity  was  silver  plated  after  assembly  to  minimize  surface  resistance.  The 
calculated  0U  was  28,000  and  the  measured  value  is  around  19,000.  The  cavity 
is  fitted  with  two  N-type  vacuum  feedthroughs  with  small  loop  antennas  for 
coupling  energy  in  and  out.  These  loops  have  thus  far  been  run  with  very 
small  coupling  coefficients,  resulting  in  =  12,000,  for  maximum  sensiti¬ 
vity.  In  future  experiments,  the  coupling  can  be  increased  to  lower  to 
about  2500,  which  will  allow  measurements  of  5  times  higher  conductivities  if 
desired. 

Sensitivity  of  the  Technique 

We  have  used  equation  (3.7)  together  with  the  v  =  constant  expressions 

for  or  and  o^  from  Table  1,  the  ( ve f £ )  values  for  Tg  =  300  K  from  Figure  1, 

DC 

and  the  cavity  parameters  given  above  to  calculate  Af/f  and  A(1  / 0)  for  a 
1-cm-diameter  centerline  plasma.  The  results  are  presented  in  Figures  6  and 
7.  Experience  to  date  indicates  that  with  the  couplings  adjusted  for  maximum 
(=  12,000),  can  detect  frequency  shifts  of  less  than  100  kHz,  which 
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FIGURE  6  FREQUENCY  SHIFT  DUE  TO  1-cm  DIAMETER  PLASMA 
IN  TM010  MODE  OF  CAVITY  AT  X  =  19.2  cm 

(T.  =  T„  =  300  K,  Rr  =  7.36  cm) 


corresponds  to  Af / f  <  6  x  10~^.  For  the  cold  plasma  conditions  of  Figure  6, 

this  corresponds  to  electron  densities  of  10*^  and  10** /cm^  for  1  amagat  of 

dry  and  humid  air,  respectively.  Likewise,  we  can  detect  Ml/Q)  changes  as 
—A  ft  o 

small  as  5  x  10”  ,  which  corresponds  to  electron  densities  of  order  10  /cm  at 

1  amagat.  Thus,  the  technique  exceeds  the  sensitivity  requirements  of  the 

task.  (It  becomes  easier  to  operate  the  experiment  with  less  sensitivity.) 
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IV  PRELIMINARY  EXPERIMENTS 


Test  Set-Up 

The  experimental  layout  for  a  set  of  experiments  designed  to  check  the 
technique  and  obtain  preliminary  data  is  shown  schematically  in  Figure  8. 

The  details  of  the  attachment  of  the  cavity  to  the  Febetron  706  are  shown  in 
Figure  9.  The  Febetron  beam  is  rather  diffuse,  and  scattering  is  accentuated 
by  the  beam's  passage  through  a  second  foil  and  hibatchi  structure,  so  that  an 
unconfined  beam  would  occupy  a  large  but  unknown  volume  of  the  cavity.  There¬ 
fore,  we  performed  most  of  these  experiments  with  a  quartz  tube,  26.5-mm  ID  by 
30-mm  OD  lining  the  cavity,  to  confine  and  define  the  area  of  the  plasma. 

This  corresponds  to  a  plasma  area  seven  times  that  of  a  1-cm-diameter  plasma, 
and  consequently  the  sensitivity  limits  in  these  experiments  corresponds  to 
electron  densities  one  seventh  of  those  quoted  above. 

The  Q  of  the  cavity  was  determined  first,  by  using  the  sweep  oscillator 
in  the  CF/AF  mode,  sweeping  ±  1  MHz  around  the  resonance  frequency.  The 
transmission  of  the  cavity  was  displayed  on  an  oscilloscope  that  had  the 
x-axis  input  driven  by  the  sweep  output.  It  was  possible  to  read  the  Af  of 
the  transmission  curve  to  about  10  kHz,  giving  the  loaded  cavity  Q  to  ~  8% 

accuracy.  The  crystal  detector  used  in  these  experiments  (HP  423B)  was  non¬ 

linear  in  the  range  of  powers  used,  and  it  was  necessary  to  correct  these  and 
other  results  for  that  effect. 

We  performed  one  experiment  that  provides  a  nice  qualitative  demonstra¬ 
tion  of  the  plasma  perturbation  technique.  We  used  an  external  ramp  function 
generator  to  frequency  modulate  the  sweep  oscillator  at  a  modulation  rate  of 

5  kHz  and  a  modulation  amplitude  of  ±  700  kHz  around  the  resonance  fre¬ 

quency.  The  transmission  history  for  several  modulation  cycles  with  an  empty 
cavity  is  shown  in  Figure  10(a).  The  transmission  curves  are  distorted  from 
the  normal  Lorentzian  shape  by  "cavity  pulling  effects,”  but  the  regularity 
and  constant  amplitude  of  the  transmission  peaks  are  clear.  Figure  10(b) 
repeats  the  transmission  history,  but  now  the  Febetron  was  fired  into  30  torr 
of  synthetic  air  in  the  cavity  at  t  ■  0.  The  first  pulse  in  the  train  has 
T  -  0  through  the  cavity,  and  the  gradual  recovery  of  peak  transmission  as  the 
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FIGURE  10  CAVITY  TRANSMISSION  FOR  30  torr  SYNTHETIC 
AIR  USING  FREQUENCY  MODULATION 
TECHNIQUE 

Upper  Trace:  Transmitted  Power 

Lower  Trace:  Frequency  Modulation  Voltage.  Shifted  So  Start  of 
Voltage  Ramp  Coincides  With  Peak  of  Unperturbed 
Cavity  Transmission 
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piasma  decays  on  a  millisecond  time  scale  is  clearly  seen.  Also,  Che  fre¬ 
quency  shifc  of  ~  30  kHz  of  the  second  pulse  and  the  smaller  shift  of  the 
third  pulse  are  clearly  detectable.  This  technique  could  be  a  rapid  way  to 
accumulate  data  if  the  results  could  be  displayed  in  a  more  accurate  format 
and  if  the  cavity  pulling  effects  could  be  eliminated  (perhaps  by  inserting  an 
isolator  in  the  line  between  the  source  and  the  cavity).  However,  it  seems 
unlikely  that  the  frequency  modulation  could  be  made  fast  enough  to  study  an 
afterglow  decay  in  full  density  air,  which  occurs  in  a  few  microseconds. 

Conductivities  at  760  and  10  torr 

The  normal  technique  for  acquiring  data  was  to  switch  the  sweep  oscil¬ 
lator  to  the  cw  mode  and  adjust  the  frequency  to  resonance,  indicated  by  a 
signal  amplitude  corresponding  to  the  peak  of  the  transmission  curve.  A 
transmission  history  at  that  frequency  was  then  acquired  using  a  normal  time 
base  sweep  triggered  by  the  Febetron  pulse.  The  source  frequency  was  then 
adjusted  to  a  number  of  +  or  -  offset  frequencies  with  the  test  repeated  at 
each  frequency.  The  sweep  oscillator  frequency  could  not  be  read  accurately 
enough  to  set  specific  frequencies,  so  instead  specific  signal  amplitudes  were 
set  (T  *  l/2Tmx,  l/^Tmax,  etc.)  and  the  frequencies  were  later  inferred  from 
the  empty  cavity  transmission  curve.  This  procedure  is  not  satisfactory 
because  of  the  uncertainty  of  the  set  frequency  and  because  the  cw  frequency 
of  this  oscillator  can  drift  over  100  kHz  in  a  moment.  Future  experiments 
will  require  a  frequency-stabilized  source. 

Representative  tracings  of  oscillographs  for  such  experiments  in  a  cavity 
opened  to  humid  laboratory  air  are  presented  in  Figure  11.  The  scope  was  set 
up  with  zero  transmission  at  the  top  of  the  figure  and  producing  nearly 

full-scale  deflection  as  in  the  top  trace.  The  signal  amplitudes  decrease  as 
the  frequency  is  tuned  off  resonance.  The  transmission  goes  to  0  at  t  *  0 
(limited  by  instrument  time  constants),  and  then  recovers  to  the  initial  value 
on  a  us  time  scale.  There  are  two  significant  features  of  these  traces. 

First  is  the  symmetry  of  the  traces  at  +  and  -  frequency  offsets,  indicating 
that  the  change  in  A(— J  is  large  compared  with  Af/f,  as  expected  at  1  amagat 
(see  Figures  6  and  7).  The  second  feature  is  the  ringing  of  the  transmission 
at  a  frequency  that  is  at  least  approximately  equal  to  the  preset  offset  fre¬ 
quency.  Understanding  and  eliminating  this  effect  will  require  further  work. 
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If  one  reads  transmissions  at  equal  time  intervals  on  each  trace  and 
crossplots  to  obtain  T  -  T(f)  at  specific  ts,  a  plot  like  Figure  12  is  ob¬ 
tained.  Neglecting  the  effects  of  ringing,  we  can  determine  three  parameters 
for  each  curve:  and*  at  *east  aPProxlmately >  ^shift^*  7,16 

first  two  parameters  offer  alternative  ways  of  determining  A(-p-](t)  from 
equations  (3.14)  and  from  the  definition  0  *  f/Af.  Determining  A(— )(t)  by  the 
two  methods  gives  two  sets  of  or(t),  as  shown  in  Figure  13.  The  two  sets  of 
results  are  in  good  agreement,  with  the  results  based  on  (T  ' /T  )  considered 
more  reliable.  Likewise,  Afs^ft(t)  can  be  used  to  determine  Cj(t),  as  shown 
in  Figure  14. 

If  one  assumes  that  these  conductivities  are  due  to  electrons,  and 

assigns  a  temperature  (or  collision  frequency)  to  them,  one  can  use  equation 

(3.7),  Table  1,  and  Figure  1  to  calculate  the  electron  densities  that  would 

cause  the  ors  and  o^s  of  Figures  13  and  14.  If  this  is  done,  the  nes 

determined  from  are  about  10  times  higher  than  those  from  or  (4  x  10^/cm^ 

9  3 

at  1  us,  compared  with  3  x  10  /cm  ).  Clearly,  something  is  amiss,  most  likely 
the  assumption  that  the  conductivities  are  due  to  electrons. 

Instead  of  assuming  a  collision  frequency,  one  can  determine  the  actual 
value  of  v  from  the  ratio  of  perturbations  (in  the  v  =  constant  approximation) 


v 


0) 

2 


(4.1) 


Results  obtained  using  the  average  A(-^)  and  Af/f  are  presented  in  Figure  15 
(division  by  p  «  amagats  does  not  affect  these  1-amagat  data). 

It  should  be  emphasized  that  the  analysis  to  this  point  has  not  required 
knowing  the  identity(s)  of  the  charged  particles.  The  conductivities  of 
Figures  13  and  14  are  correct,  as  are  the  vs  of  Figure  15.  If  we  assume  that 
they  are  electrons,  we  can  calculate  ng(t)  as  shown  in  Figure  16.  The  results 
from  the  three  data  sources  are  now  in  good  agreement,  as  indeed  they  must  be 
since  the  experimental  v  was  used  in  the  calculations. 
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FIGURE  13  REAL  CONDUCTIVITY  HISTORIES  FOR  1  amagat 
HUMID  AIR  (Rp  =  1.3  cm) 
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FIGURE  14  IMAGINARY  CONDUCTIVITY  HISTORY  FOR  1 
HUMID  AIR  (R„  =  1.3  cm) 
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FIGURE  16  DEDUCED  ELECTRON  DENSITY  HISTORIES  FOR 
1  amagat  HUMID  AIR 


It  is  unlikely  that  electrons  are  the  charge  carrier  in  these  measure¬ 
ments  for  two  reasons.  First,  a  comparison  of  the  vs  from  Figure  15  with  the 
calculated  values  shown  in  Figure  1  shows  that  the  measured  values  are  more 
than  10  times  smaller  than  the  calculated  values  for  humid  air.  Second,  the 
attachment  rate  of  electrons  to  0^  in  air  at  1  amagat  is  very  fast, 


-dneN2  =  -g'ne 


(4.2) 


7  7  —1 

where  8'  =  8  x  10  p  s  for  room  temperature  electrons.  This  attachment 
rate  should  reduce  the  electron  density  below  104/cm^  by  300  ns,  whereas 

Q  O 

Figure  16  indicates  ng  ■  10  /cm  at  2  us.  The  microwave  fields  in  the  cavity 
might  heat  the  electrons  and  slow  down  attachment  somewhat,  but  not  enough  to 
account  for  such  differences. 

The  other  identification  of  charged  particles  is  ions,  now  viewed  as 
positive /negative  ion  pairs.  Ion  conductivities  are  usually  neglected  because 
when  ne  =  N^, 


(_E£.j  =  -L  >  5.8  x  104  (for  O?")  (4.3) 

vui  ,  m  * 

pi  e 


(where  m  is  the  mass)  so  that  electrons  dominate  the  imaginary  conductivity, 
and  the  mobility  ratios  for  electrons  and  ions  are 

ye  4 

~  »  104  (4.4) 

M1 

so  that  again  electrons  dominate  the  real  conductivity.  This  neglect  of  ions 
is  no  longer  valid  when  n^  »  104  ng. 

If  we  assign  the  conductivity  to  ions,  then  with  an  assumption  about  m^, 
we  can  calculate  n^  since  we  have  already  determined  v  experimentally.  Let  us 
assume,  for  example,  that  the  dominant  negative  ion  is  0^”  and  the  dominant 
positive  ion  is  H.jO+,H20.  Then  the  effective  mass  (in  atomic  units)  is 
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m 


al°2 


eff  m,  + 


1  +  “2 


64  x  37 
64  +  37 


23 


(4.5) 


The  plasm  frequency,  u)p  »  2 it  f  is  then  related  to  the  Ion  density  by 


n.e 


“pl 


“  J  ^ 

- ) 

meff£o 


1/2  7  ,  ,n4  1/2 

=  7.6  x  10  n^ 


(4,6) 


We  can  then  use  the  v  =  constant  expression  for  either  ar  or  Oj  from  Table  1 

to  calculate  u)p  and  consequently  n^  from  the  results  of  Figure  13  or  14.  The 

resulting  ion  density  history  is  presented  in  Figure  17.  The  ion  densities 

are  4  x  10^*  times  larger  than  the  electron  densities  of  Figure  16;  this  is 

just  the  ratio  m  ff/ta 
err  €• 

The  collision  frequency,  v,  is  related  to  the  collisional  cross  section, 
ac,  and  the  drift  velocity,  vD,  by  the  relationship 

v  «=  nocvD  .  (4.7) 


If  we  approximate  vD  by  the  average  Ion  velocity 

vD  =  v  (8kT/nmeff )1^2  ,  (4.8) 

which  is  v  a  5  x  10^  cm/s,  we  find  that  ac  s  1.3  x  10  ^  cm^.  This  is  only 
somewhat  larger  than  a  typical  ion  collision  cross  section  of  oc  *  5  x 
10'1-*  cm^  [McD64],  perhaps  reflecting  the  fact  that  vQ  is  larger  than  v  .  The 
collision  frequency  can  also  be  related  to  the  ion  mobility  K  by 


K 


2  2 
m  ,,e  u>  +  v 
eff  o 


(4.9) 


2 

This  yields  K  =  2.0  cm  /V-s,  which  is  in  excellent  agreement  with  typical 
values  for  ions  at  low  electric  fields  {McD64}.  Both  these  calculations 
support  the  assignment  of  the  measured  conductivity  to  ions  for  this  1  amagat 
humid  air  case. 
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FIGURE  17  DEDUCED  ION  DENSITY  HISTORY  FOR  1  amagat 
HUMID  AIR 

(Assume  =  23  amu) 
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For  10-torr  synthetic  air,  the  situation  seems  much  more  straightfor¬ 
ward.  Representative  scope  traces  of  cavity  transmission  at  fQ  and  three  ±Afs 
are  shown  in  Figure  18;  there  are  three  marked  differences  from  the  760  torr 
data  of  Figure  12.  First,  the  time  scale  for  plasma  decay  is  about  250  times 
as  long.  Second,  the  +  and  -  Af  curves  are  now  unsymuKtric  because  the 
frequency  shift  is  larger  than  the  broadening  due  to  a(-^),  as  expected  from 
Figures  6  and  7.  Third,  the  superimposed  ringing  is  no  longer  present. 
Crossplotted  signal  histories  at  250  us  inte^vqls  are  presen|ed  in 
Figure  19.  Again,  we  can  determine  or  from  (t)  or  from  (^)(t)  »  with  the 
results  shown  in  Figure  20.  Again,  the  two  techniques  give  good  agreement. 

If  we  now  calculate  v  from  equation  (3.22)  and  normalize  the  results  to  1 
amagat,  we  obtain  the  results  also  shown  in  Figure  15.  These  collision 
frequencies  are  consistent  with  the  dry  air  results  of  Figure  1  at 
approximately  0.15  eV  (=  1700  K).  The  electrons  are  probably  in  equilibrium 
with  the  vibrational  mode  of  N2  at  this  temperature.  It  is  less  likely  that 
the  elevated  Te  is  due  to  the  microwave  field  heating  the  electrons,  but  that 
effect  should  be  checked  by  varying  the  source  power. 

The  electron  density  histories  deduced  from  ar,  a4,  and  v  are  presented 
in  Figure  21.  Presumably,  the  decay  is  due  to  attachment,  which  finally  takes 
over  when  the  densities  drop  too  low  for  recombination  to  dominate.  Such  a 
decay  would  be  exponential  (linear  on  Figure  21).  The  best  linear  fit  to  the 
data  gives  8*  =  2  x  10  p  s  ,  compared  with  a  value  for  room  temperature 
electrons  of  S'  5  8  x  10  p  s  .  A  higher  electron  temperature  of  “  0.15  eV 
deduced  from  the  v's  of  Figure  15  would  be  expected  to  increase  the  attachment 
rate  by  a  factor  of  4  [CPB62],  in  contrast  with  the  reduction  observed  in  the 
present  results.  The  reduction  and  the  nonexponential  electron  decay  may  be 
due  to  detachment  of  electrons  from  negative  ions  by  collision  with  excited 
species,  which  could  include  Nj  (vib),  ^(A),  02(*A),  and  02(*Z). 

Pressure  Scans 

The  results  reported  above  involved  15  to  20  tests  in  each  case,  with  the 
attendant  amount  of  data  reduction.  It  was  not  appropriate  to  repeat  this 
lengthy  procedure  for  preliminary  tests  at  many  pressures,  but  we  did  make 
several  measurements  at  fQ  only  over  a  wide  pressure  range.  These  results  can 
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FIGURE  20  REAL  CONDUCTIVITY  HISTORIES  FOR  10  torr 
SYNTHETIC  AIR 
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be  used  to  determine  <Jr  by  the  (T  ' /T  )  method  if  we  neglect  the  fact  that 
part  of  the  reduction  of  transmission  at  f  ■  fQ  results  from  a  frequency  shift 
rather  than  a  A(l/Q)  broadening.  This  approximation  will  cause  the  most 
severe  errors  In  low  pressure  and  early  time  results. 

Transmission  histories  for  two  series  of  tests  that  compare  the  results 
of  laboratory  and  synthetic  air  are  presented  in  Figure  22.  Except  for  the 
p  -  760  torr  results,  the  conductivity  decays  are  always  significantly  longer 
in  synthetic  air  than  in  laboratory  air.  It  is  known  that  HjO  Is  a  much  more 
effective  third  body  in  attachment  of  electrons  to  O2 ,  and  the  results  of 
Figure  22  can  be  explained  if  the  partial  pressure  of  F^O  stayed  constant  as 
the  total  pressure  was  reduced.  We  also  find  that  the  decay  times  in  synthe¬ 
tic  air  rapidly  decrease  with  impurity  level,  causing  considerable  scatter  in 
these  preliminary  results. 

The  conductivity  decays  get  progressively  longer  with  decreasing  pres¬ 
sure,  until  at  3  torr  in  synthetic  air  and  1  torr  in  laboratory  air,  the 
character  of  the  decay  changes  markedly.  Me  assume  that  these  are  the  pres¬ 
sures  where  diffusion  of  electrons  to  the  walls  becomes  faster  than  attach¬ 
ment.  These  rather  low  pressures  then  mark  the  lower  limits  of  these  measure¬ 
ments  in  air. 

We  pointed  out  above  that  we  can  determine  ar(t)  from  these  traces  if  we 
neglect  the  contribution  of  frequency  shifts  to  reduction  of  transmission. 

This  neglect  would  tend  to  make  the  or  values  low  at  the  beginning  of  each 
trace,  particularly  at  lower  pressures.  Results  for  two  series  of  tests  in 
synthetic  air,  one  with  a  quartz  tube  of  26.5-mm  ID  lining  the  cavity  and  the 
second  using  a  20-mm-ID  tube,  are  presented  in  Figures  23  and  24.  The  purpose 
of  the  two  series  was  to  look  for  a  possible  wall  effect  on  the  measure¬ 
ments.  The  two  sets  of  results  are  similar  but  not  identical,  and  the 
accuracy  of  these  preliminary  data  is  not  good  enough  to  warrant  a  firm  con¬ 
clusion.  A  third  series  of  tests  with  no  tube  at  all  showed  quite  similar 
conductivity  decay  times,  but  since  we  do  not  know  the  plasma  area  in  that 
case,  we  cannot  calculate  the  ors. 
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Laboratory  Air 


FIGURE  22  CAVITY  TRANSMISSION  HISTORIES  FOR  SYNTHETIC  AND 
LABORATORY  AIR 

(f  =  f0  =  1.537  GHz,  plasma  diameter  =  2.0  cm,  cavity  band¬ 
width  "•  160  kHz) 
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FIGURE  22  CAVITY  TRANSMISSION  HISTORIES  (Concluded) 
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V  CONCLUSIONS  AND  PLANS 


The  work  reported  here  indicates  that  the  microwave  cavity  perturbation 
technique  will  be  a  powerful  tool  for  diagnosing  the  afterglow  channel  of  high 
energy  electron  beams.  The  most  sensitive  measurement  is  of  the  real  conduc¬ 
tivity  (which  is  the  D.C.  conductivity  at  high  pressures).  We  have  demon¬ 
strated  detectivity  limits  of  ar  >  10^  s-*  for  a  2.65-cm-diameter  plasmas;  the 
sensitivity  varies  inversely  with  plasma  diameter,  so  that  the  threshold  for  a 
1-cm-diameter  plasma  is  or  >  10^  s-* .  This  is,  in  fact,  somewhat  more 
sensitive  than  desired  for  diagnosing  ETA  or  ATA  experiments,  but  the  sensiti¬ 
vity  can  easily  be  reduced  (most  readily  by  increasing  the  coupling  of  the 
input  and  output  loops).  Thus  far,  our  sensitivity  limit  for  the  imaginary 
conductivity  has  been  about  the  same  as  that  of  the  real  conductivity,  or 
Oj  >  5  x  10^  s"1  (again  for  a  2.65-cm-diameter  plasma). 

Either  the  real  or  imaginary  conductivities  can  be  used  to  deduce 
electron  densities  if  a  collision  frequency  is  assumed.  Alternately,  the 
ratio  of  conductivities  gives  the  collision  frequency  directly,  and  the 
electron  densities  can  then  be  calculated.  Furthermore,  comparing  the 
measured  collision  frequencies  with  calculated  values  allows  one  to  deduce 
electron  temperatures.  Since  electron  energies  are  often  equilibrated  with 
vibrational  energies,  this  indicates  vibrational  temperatures  as  well.  Thus, 
the  technique  offers  determination  of  a  wide  range  of  information. 

When  attachment  dominates  over  recombination  in  the  early  afterglow, 
negative  ion  densities  can  become  several  orders  of  magnitude  larger  than 
electron  densities.  Then,  the  conductivity  is  due  to  ions  rather  than 
electrons.  In  our  measurements  to  date,  this  effect  manifested  itself  when 
collision  frequencies  were  inconsistent  with  values  expected  for  electrons. 

The  conductivity  measurement  remains  valid  even  in  this  case.  Furthermore, 
application  of  pulsed  fields  in  the  cavity  (microwave  or  D.C.)  to  detach  bound 
electrons  from  ions  might  simulate  the  fields  induced  by  a  high  current  elec¬ 
tron  beam,  and  furnish  information  on  the  change  in  residual  conductivity  that 
occurs  in  the  nose  of  a  second  or  subsequent  pulse. 
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There  are  many  refinements  to  be  made  to  the  system  before  taking  final 
data.  First,  a  much  more  stable  microwave  source  is  needed  so  that  recordings 
can  be  taken  at  fixed  and  known  frequencies  (this  source  is  currently  being 
procured).  A  high  purity  vacuum  system  is  needed  to  take  data  with  synthetic 
air  because  the  preliminary  results  showed  that  the  conductivity  was  sensitive 
to  impurities  (presumably  water  vapor).  It  would  be  helpful  to  add  a  magnetic 
field  to  the  cavity  to  confine  the  Febetron  electron  beam  to  a  more  nearly 
cylindrical  geometry.  With  or  without  the  field,  the  beam  axial  current 
profile  should  be  mapped  with  a  Faraday  cup  probe  and/or  open  shutter 
photography.  Finally,  the  effect  of  microwave  power  level  on  the  afterflow 
decay  should  be  investigated. 

These  refinements  will  lead  to  measurements  using  our  Febetron  706 
electron  beam  source  that  should  provide  definitive  conductivity  histories  for 
room  temperature  air  as  a  function  of  pressure.  The  experiment  is  also 
designed  to  fit  directly  onto  the  ETA  facility  for  similar  measurements  on 
that  beam.  This  will  be  particularly  interesting  when  ETA  is  operated  in  the 
burst  pulse  mode.  Future  work  will  also  include  extending  the  measurements  to 
high  temperature  (shock-heated)  air  to  investigate  the  effects  of  temperature 
on  conductivity  decay. 
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